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ABSTRACT 

We estimate and remove the contamination of weak gravitational lensing measurements by 
the intrinsic alignment of close pairs of galaxies. We do this by investigating both the aper- 
ture mass B mode statistic, and the shear correlations of close and distant pairs of galaxies. 
These can be used to quantify non-lensing effects in weak lensing surveys. We re-analyse 
the COMBO-17 survey, and study published results from the Red-sequence Cluster Survey 
and the VIRMOS-DESCART survey, concluding that the intrinsic alignment effect is at the 
lower end of the range of theoretical predictions. We also revisit this theoretical issue, and 
show that misalignment of baryon and halo angular momenta may be an important effect 
which can reduce the intrinsic ellipticity correlations estimated from numerical simulations 
| to the level that we and the SuperCOSMOS survey observe. We re-examine the cosmological 

. parameter estimation from the COMBO-17 survey, using the shear correlation function, and 

now marginalising over the Hubble constant. Assuming no evolution in galaxy clustering, and 
marginalising over the intrinsic alignment signal, we find the mass clustering amplitude is 
reduced by 0.03 to <j$(£l m /0.27) ' 6 — 0.71 ± 0.11, where f2 m is the matter density param- 
eter. We consider the forthcoming SuperNova/Acceleration Probe wide weak lensing survey 
(SNAP), and the Canada-France-Hawaii Telescope Legacy wide synoptic survey, and expect 
them to be contaminated on scales > 1 arcmin by intrinsic alignments at the level of ~ 1% 
and ~ 2% respectively. Division of the SNAP survey for lensing tomography significantly 
increases the contamination in the lowest redshift bin to ~ 7% and possibly higher. Removal 
of the intrinsic alignment effect by the downweighting of nearby galaxy pairs will therefore 
be vital for SNAP. 

b : 

Key words: cosmology: observations - gravitational lensing - large scale structure, galaxies: 
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1 INTRODUCTION 

The detection of weak gravitational lensing by large scale structure 
is a direct way to measure the total matter distribution in the Uni- 
verse, demanding no assumptions for how luminous matter traces 
the dominant, largely unknown, dark matter component. Detected 
by several groups, weak gravitational lensing is now a well es- 
tablished technique, used successfully to set joint constraints on 
the matter density parameter $l m and the amplitude of the mat- 
ter power spectrum, as, (R.Maoli et al. 2001; Rhodes, Refregier 
& Groth 2001; Van Waerbeke et al. 2001; Hoekstra, Yee & Glad- 
ders 2002; Bacon et al. 2003; Jarvis et al. 2003; Brown et al. 2003; 
Hamana et al. 2002), to measure the bias parameter b (Hoekstra, 
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Yee & Gladders 2001; Pen et al. 2003), and has recently been used 
to directly extract the 3D non-linear matter power spectrum Ps(k) 
(Tegmark & Zaldarriaga 2002; Pen et al. 2003). Combined with 
cosmic microwave background observations, weak lensing can pro- 
vide strong constraints for as and fl m as the degeneracies in each 
measurement are almost orthogonal in the ag, — Q m plane (Brown 
et al. 2003; Contaldi, Hoekstra & Lewis 2003). 

Unlike many other tests of cosmology, weak lensing surveys 
with photometric redshift information can tightly constrain cosmo- 
logical quintessence models and the equation of state parameter 
w, which will be key to our understanding of dark energy, (Heav- 
ens 2003; Refregier et al. 2003; Benabed & Van Waerbeke 2003; 
Jain & Taylor 2003). With the increased image resolution avail- 
able from multi-colour space-based lensing surveys it will also be 
possible to construct high-resolution projected dark matter maps, 
and 3D dark matter maps of mass concentrations > 1 x 1O 13 M0 
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(Taylor 2001; Hu & Keeton 2002; Bacon & Taylor 2003; Massey 
et al. 2003). With future deeper and wider multi-colour surveys, 
for example the Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS) (www.cfht.hawaii.edu/Science/CFHTLS) and the space- 
based SuperNova/ Acceleration Probe (SNAP) (snap.lbl.gov), weak 
gravitational lensing will soon reach its 'era of high precision cos- 
mology', provided it can get a good handle on the many causes of 
systematic errors that arise when trying to detect this minute weak 
lensing signal. 

Sources of systematic errors in weak lensing analysis arise 
from the shearing and smearing of galaxy images caused by the 
atmosphere, telescope optics, and detectors (Kaiser, Squires & 
Broadhurst 1995; hereafter KSB; Luppino & Kaiser 1997). With 
excellent seeing observing conditions, or space-based data, com- 
bined with instruments and detectors that are designed with weak 
lensing detection in mind, these effects can be minimised and cor- 
rected for (see for example Rhodes et al. 2003 and Bacon, Refregier 
& Ellis 2000). Aside from these observational effects there is a po- 
tentially significant error arising from a key assumption made for 
all weak lensing studies, that galaxy ellipticities are randomly ori- 
ented on the sky. Gravitational interactions during galaxy forma- 
tion could however produce intrinsic shape correlations between 
nearby galaxies, mimicking to an extent weak lensing shear cor- 
relations. As the new generation of wide-field deep weak lensing 
surveys beat down their observational sources of systematics, it is 
this additional source of intrinsic ellipticity correlations that could 
limit the accuracy of cosmological parameter estimation from weak 
lensing studies. The extent to which this is true will be aided by the 
study in this paper. 

Observational evidence for the existence of intrinsic galaxy 
alignments comes from the detection of galaxy ellipticity corre- 
lations in low redshift surveys where weak lensing shear corre- 
lations are negligible, for example in the SuperCOSMOS survey 
(Brown et al. 2002), and in the Tully catalogue (Lee, Pen & Sel- 
jak 2002). Theoretically, the intrinsic alignment of nearby galaxies 
has been investigated through numerical simulations and analyti- 
cal techniques. These have provided estimates of the order 10% 
contamination to weak lensing measurements from surveys with 
median redshift z m ~ 1.0 (Heavens, Refregier & Heymans 2000, 
hereafter HRH, Croft & Metzler 2000; Catelan, Kamionkowski & 
Blandford 2001; Crittenden et al. 2001, hereafter CNPT; Lee & 
Pen 2001; Hui & Zhang 2002; Porciani, Dekel & Hoffman 2002; 
Jing 2002, hereafter Jing; Mackey, White & Kamionkowski 2002). 
Whilst there is broad agreement between these studies on the effect 
for weak lensing measurements, the finer details can differ by up 
to an order of magnitude or more, with the numerical simulations 
generally predicting a higher level of contamination than the semi- 
analytic studies. In this paper, we re-examine this issue, and show 
that a combination of misalignment of the baryon and halo angular 
momentum, as determined by van den Bosch et al. (2002), and the 
finite thickness of disk galaxies, modifies the predictions of HRH 
bringing them into good agreement with the semi-analytic model 
of CNPT. 

With redshift information, it has been shown that the intrinsic 
signal can be suppressed in weak lensing analysis by downweight- 
ing galaxy pairs which are physically close (Heymans & Heavens 
2003; King & Schneider 2002). This can be done optimally with- 
out significantly increasing the shot noise in the final weak lensing 
analysis, and for this it is helpful (although not necessary) to have 
a good estimate of the level of contamination. Obtaining an ob- 
servationally constrained estimate is therefore one of the purposes 
of this paper, using three weak lensing surveys: a re-analysis of 



COMBO- 17 (Brown et al. 2003), and the published results of the 
Red-sequence Cluster Survey (RCS; Hoekstra et al. 2002) and the 
VIRMOS-DESCART survey (Van Waerbeke et al. 2002). 

This paper is organised as follows. In Section 2 we review 
three different theoretical models for the intrinsic ellipticity corre- 
lations between nearby galaxies, taken from HRH, Jing and CNPT. 
We present a modification to the HRH analysis that we apply to 
numerical simulations in Section 2.2, finding excellent agreement 
with the observed intrinsic alignment signal from the SuperCOS- 
MOS survey. Using the three different intrinsic alignment models: 
HRH, Jing, CNPT, and our modified HRH model which we shall 
call HRH* hereafter, we then determine intrinsic alignment con- 
tributions to the aperture mass B mode statistic M± (Schneider 
et al. 1998), for the RCS and the VIRMOS-DESCART survey. Us- 
ing published measurements of M± as upper limits for the intrinsic 
galaxy alignment contribution, we show in Section 3 that, assuming 
there is no evolution in galaxy clustering, we can reject the intrin- 
sic alignment models of Jing and HRH. In Section 4, we obser- 
vationally constrain and remove the contribution to COMBO- 17's 
weak lensing measurements by intrinsic galaxy alignments. This is 
made possible due to the highly accurate photometric information 
of the COMB O-l 7 survey, through the application of an optimal 
intrinsic alignment contamination removal method as described in 
Heymans & Heavens (2003). In Section 5, we investigate the ef- 
fect that galaxy clustering evolution would have on our results. We 
look at the implications for weak lensing analysis in Section 6, con- 
straining ag and Q. m with the shear correlation function statistic 
from COMBO- 17, where we now include marginalisation over our 
measured intrinsic alignment signal. We also determine estimates 
of the contamination of shear correlation measurements from the 
CFHTLS and SNAP surveys, summarising and concluding in Sec- 
tion 7. 



2 INTRINSIC ELLIPTICITY CORRELATIONS 

The ellipticity of a galaxy approximated as an ellipse with axial 
ratio q, at position angle <f> measured counter-clockwise from the x 
axis, can be defined as 

( ei \ g 2 -l ( cos2</) \ 

\ £2 J q 2 + l\ sm2(j) J ' W 

In the weak lensing regime the observed (complex) ellipticity e is 
related to the source ellipticity e s and the complex shear 7 by e ~ 
e a + 2~/. With many source galaxies the shear can be estimated from 

7- \{e) = i(ei + ie 2 ). (2) 

See Rhodes, Refregier & Groth (2000) for an exact relationship. 
A measurement of the shear correlation function (77*) can be 
directly related to the matter power spectrum (see Bartelmann & 
Schneider 2001 and references therein), but to estimate it observa- 
tionally includes an uncertain contribution from intrinsic correla- 
tions between the ellipticities of nearby galaxy pairs I a b = e" a el* 
if, when averaging over many galaxy pairs, (I a b) 7^ 0. Hence the 
expectation of the shear correlation function is given by 

£[77* ■■e] = lYl M X K( X + °) - ^) ■ ( 3 ) 

ab 

Several groups have found theoretical support for non-zero intrinsic 
correlations, and the existing methods used to determine its ampli- 
tude can be broadly split into two approaches. The first employs 
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numerical simulations with varying simplifying assumptions, relat- 
ing galaxy shape to dark matter halo properties. The second uses 
analytic techniques, relating ellipticity correlations to initial linear 
tidal shear field correlations. These methods approximate I ab as 
a theoretically determined three-dimensional ellipticity correlation 
function rj(r) which depends only on the comoving distance be- 
tween the galaxy pair r = |r|: 



r)(r) = (e(x)e*(x + r)). 



(4) 



We project rj(r) into two dimensions using a modified Limber 
equation, in order to compare it with observed angular elliptic- 
ity correlations. In the absence of weak gravitational lensing, for 
example with the low redshift SuperCOSMOS data (Brown et al. 
2002), measured angular ellipticity correlations can be directly 
compared to those predicted from theoretical models: 



Clip) = 



J dz a dz b <t> z (z a )4> z (z b ) [1 + igg{r ab )]r](r ab ) 

j dz a dz b (f) z (Za)4>z(z b ) [1 + £gg(r a b)] 



(5) 



where <f> z (z a ) is either a known redshift distribution, for example 
the COMBO- 17 photometric redshift distribution shown in Fig. Al 
of the Appendix, or can be approximated by 



z (z) oc z 2 exp 



(6) 



for a survey with median redshift z m w 1.4?o. Most studies have 
found T](r a b) to be significantly non-zero only for galaxies closer 
than a few tens of Mpc, where the comoving separation, r ab is given 
accurately enough in a flat universe by 

2 



r 2 ■ 

•ab 



(W a — Wb) + 



w a + w b 



(7) 



where w is the comoving radial geodesic distance. The signal de- 
pends on the galaxy two-point correlation function, £ gg (r), since it 
determines how many galaxy pairs observed at a given angular sep- 
aration are physically close together. We take the correlation func- 
tion to be non-evolving (Cole et al. 1994; Giavalisco et al. 1998; 
Adelberger et al. 1998): 



(8) 



with 7 = —1.8 and ro = 5/i _1 Mpc. In Section 5 we will also 
consider the effect of redshift dependence in the galaxy correla- 
tion function, as claimed by Le Fevre et al. (1996); Postman et al. 
(1998); Baugh et al. (1999); Carlberg et al. (2000); Firth et al. 
(2002); Wilson (2003). 

In this Section we will focus on three studies: HRH and Jing, 
which are based on numerical simulations, and CNPT, where the 
strength of ellipticity correlations are analytically computed in lin- 
ear theory. Fig. 1 compares three published functional fits, de- 
fined in Section 2. 1 , for the three-dimensional ellipticity correlation 
functions, r)(r), showing that whilst there is some rough agreement 
in the results from these studies, the numerical simulations in gen- 
eral, predict a stronger intrinsic alignment effect compared to the 
analytical technique. 

2.1 Intrinsic Alignment Models 

Assuming that luminous matter forms galaxies in all dark matter 
halos above a minimum mass limit, Mh, it is possible to acquire a 
three-dimensional catalogue of galaxy shapes from a large N-body 
dark matter simulation. HRH have modelled the shape of each lu- 
minous galaxy as an infinitesimally thin disc placed perpendicular 
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Figure 1. Intrinsic ellipticity correlation models. Comparison of differ- 
ent functional fits to three intrinsic alignment studies, Jing (solid), HRH 
(dashed) and CNPT (dot dash). The modified HRH model that is measured 
and defined in Section 2.2 is also plotted (dotted). 



to its angular momentum vector, which is assumed to be perfectly 
aligned with that of its parent dark matter halo. This resulted in 
ellipticity correlations that were fitted with an exponential 



J7hrh(V) = 0.012 exp - 



1.5ft" 1 Mpc 



(9) 



Using a higher resolution N-body simulation, Jing assumed that the 
galaxy ellipticity was equal to that of its parent dark matter halo, re- 
sulting in ellipticity correlations with a best fitting functional form 
defined for (erei), 



3.6 x 10~ 2 



??jin g (r) « 2(eiei) = 2 



( Mi V 

^ 10 10 h- 1 Mpc J 



r°- 4 (7.5 17 +r 1 - 7 ) 



(10) 



where in this paper we will use Mh = 6.9 x 10 11 M Q , which cor- 
responds to the minimum mass limit used by HRH. HRH found a 
similar shape for the ellipticity correlations of dark matter haloes. 
CNPT related galaxy ellipticity correlations to the initial correla- 
tions in the tidal shear field. Assuming that the density field has a 
Gaussian distribution, CNPT found that the ellipticity correlation is 
well approximated at large r > 3ft -1 Mpc by 



J7cnpt(V) = 



a 2 (3 2 e\r) 
84 e 2 (0)' 



(11) 



where a quantifies the uncertainty in the correlation between the 
shear field and a galaxy's moment of inertia, (3 accounts for the non- 
zero thickness of galaxy disks where e ga i = /3e t hindisc, and e(r) is 
a density correlation function. For the purpose of this paper, we use 
an analytic solution to equation (1 1) for the simple model proposed 
in CNPT, who took a = 0.55, inferred from numerical simulation 
results (Heavens, Refregier & Heymans 2000), /3 = 0.73, from the 
mean of the measured ellipticity distribution of the APM survey 
(Lambas, Maddox & Loveday 1992) and from the observed distri- 
bution of ellipticities measured for lensing studies (Ebbels 1998), 
and e(r) oc 1/r, where we also follow CNPT by including top-hat 
smoothing on l/i _1 Mpc scales. 

2.2 Intrinsic alignments of galaxies: revisiting theoretical 
expectations 

Recent results from numerical simulations of structure formation 
by van den Bosch et al. (2002), have shown that the angular mo- 
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Figure 2. Distribution of the misalignment angle between the total angular 
momentum vectors of the dark matter and gas, as found by van den Bosch 
et al. (2002). Over-plotted is the best fit Gaussian model with [m = 19.5° 
and cr = 19.3. 
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Figure 3. 3D ellipticity correlation function estimated from the ACDM 
Virgo simulation, assuming that the gas disc and dark matter halo are mis- 
aligned. Over-plotted is the HRH model ijhrh ( r ) (dotted) and the best fit 
HRH modified model r/ H R, H * ( r ) (dashed). 



mentum of non-radiative gas and dark matter are actually poorly 
aligned with a median misalignment angle of ~ 30°. Observational 
support for some misalignment, albeit smaller, comes from strong 
lensing considerations in the CASTLES survey (Keeton, Kochanek 
& Falco 1998). With these results we are able to advance the analy- 
sis of HRH, where instead of assuming perfect alignment between 
the angular momentum of the halo and baryons, we include a prob- 
ability distribution for the misalignment. Figure 2 shows the prob- 
ability distribution of the polar misalignment angle 9 as found by 
van den Bosch et al. (2002) and the best fit Gaussian model, van 
den Bosch et al. (2002) note that higher resolution simulations are 
required to reduce the possible impact of discreteness effects on 
this probability distribution, but they do conclude that there is a 
true misalignment that is significantly different from zero. 

In light of these results we have reanalysed the N-body simu- 
lation developed by the Virgo Consortium (Jenkins et al. 1998), for 
a ACDM cosmological model, (Q m = 0.3, Qa = 0.7, a 8 = 0.9, 
r = 0.2), approximating the probability distribution for the mis- 
alignment polar angle as found by van den Bosch et al. (2002), 
as a Gaussian with mean fi — 20° and width a — 20°, which 
is truncated at zero misalignment. We assign a random misalign- 
ment azimuthal angle around the original halo angular momentum 
vector. The galaxy ellipticity is then determined from the angular 
momentum of the baryons, by, 



HL z )\=p 



1 + Li 



(12) 



where we now account for galaxy thickness, following CNPT by 
taking j3 = 0.73. Fig. 3 compares the resulting 3D ellipticity cor- 
relation function r){r), with the HRH model from equation (9), 
(dashed), showing that with the inclusion of baryon halo angular 
momentum misalignment, the amplitude of the galaxy alignments 
predicted from numerical simulations is significantly reduced. We 
choose to fit a simple function 



r){r) = 



A 



1 + (r/B) 2 



(13) 



We fix B — lh^ 1 Mpc, which is similar to the CNPT choice of a 
smoothing length of lft _1 Mpc, and calculate the maximum likeli- 
hood value for the amplitude A. A small positive value is preferred, 
although as shown dot-dashed in Fig. 10, it should be noted that the 
reduced signal from the misalignment means the data are actually 



consistent with zero. The best fit shown dashed in Fig. 3 is (r in 
h- 1 Mpc) 



0.0011 



(14) 



l +r 2 ' 

which is consistent with the CNPT analytical model. Note that the 
new model fit lies above the old HRH model at large scales. This 
should not be taken to imply that misalignment increases the corre- 
lations, rather, it reflects the fact that the old exponential fit of HRH 
underestimated the (noisy) correlations on large scales. 

Comparing our HRH* model with the ellipticity variance mea- 
sured in the low redshift SuperCOSMOS survey (Brown et al. 
2002), we find excellent agreement, especially on scales > 35 ar- 
cmin, see Fig. 4. We have calculated the HRH* angular elliptic- 
ity correlation function from equation (5) with a median redshift 
z m — 0.1, and related it to the ellipticity variance statistic follow- 
ing Bacon et al. (2003): 



<7*(0), 



20 F 

e 2 



H 



de' e'Ci(e'), 



(15) 



where the factor ^ is a good approximation to use in order to 
scale the ellipticity variance measured in circular apertures of ra- 
dius 9 to the ellipticity variance as measured in SuperCOSMOS, in 
square cells of length 9 (Bacon, Refregier & Ellis 2000). We cal- 
culate the likelihood of the HRH* model from the SuperCOSMOS 
results on scales 9 > 20 arcmin, allowing the amplitude to vary, see 
Fig. 10. Below this angular scale the SuperCOSMOS results are ex- 
tremely noisy and could suffer contamination from residual point 
spread function anisotropy distortions, (see Brown et al. (2002) for 
details). We therefore omit these points from the likelihood fitting 
procedure, finding the best fit amplitude A = 0.0009 ± 0.0005, 
(95% confidence limits). 



3 UPPER LIMITS ON THE INTRINSIC ALIGNMENT 
SIGNAL: APERTURE MASS B MODES 

A good diagnostic for determining the level of systematic errors 
present in weak lensing measurements, is to decompose the shear 
correlation signal into E and B modes. This was first proposed by 
Crittenden et al. (2002), and is now a standard statistical test for 
the presence of non-lensing contributions to weak lensing measure- 
ments. Weak gravitational lensing produces curl-free distortions 



© 0000 RAS, MNRAS 000, 000-000 



Intrinsic Galaxy Alignments: estimation and removal 5 




60 

S(arcmins) 



100 



Figure 4. Comparison of the measured ellipticity variance in the Super- 
COSMOS survey, see Brown et al. (2002) for details, compared to the el- 
lipticity variance expected at this redshift from the modified HRH model, 
HRH*. 



(E-type), and contributes only to the B-type distortions at small 
angular scales, 6 < 1 arcmin, due to source redshift clustering 
(Schneider, Van Waerbeke & Mellier 2002). A significant detection 
of a B-type signal in weak lensing surveys is therefore an indica- 
tion that ellipticity correlations exist either from residual systemat- 
ics within the data and/or from intrinsic galaxy alignments. 

The decomposition has previously been carried out by de- 
termining the B mode shear power spectrum, Cf 3 (Brown et al. 
2003; Pen, Van Waerbeke & Mellier 2002) or by using the aperture 
mass statistic, M ap (Hoekstra, Yee & Gladders 2002; Van Waer- 
beke et al. 2002; Hamana et al. 2002; Jarvis et al. 2003), where 
the aperture mass statistic is defined to be the aperture weighted di- 
mensionless surface mass density of a lens, k (Kaiser et al. 1994), 
such that 



M ap (0) = / d'#U{&)K(d) , (16) 
Jo 

where the simplest radial filter function is given by, 



(Schneider et al. 1998). With this choice of filter function the aper- 
ture mass Map is related to the power spectrum by 



(Mi P )(8) 



n9 4 



dk 

F 



P K (k) [Ji{k0)f 



(18) 



where ^4,(k6) is the fourth-order Bessel function of the first kind, 
and P re (fc) is the convergence power spectrum at wave number k, 
related to the non-linear mass power spectrum Ps by 



P*{k) = 



9H^ 

4c 4 



dwCr^Ps 



a 2 (w) \ }k (w) 



k 



(19) 



a(w) is the dimensionless scale factor, Hq is the Hubble param- 
eter and tt m the matter density parameter. The second argument 
of Ps allows for time-evolution of the power spectrum. g(w) is a 
weighting function locating the lensed sources, 



g(w) = 



j / ,/ is fa (w - w) 



(20) 



<f)(w(z))dw is the observed number of galaxies in dw, and wh is 
the horizon distance (Schneider et al. 1998). 

In this paper we will focus on the aperture mass decomposi- 
tion, which can be directly calculated from angular ellipticity cor- 
relation functions as follows: 



B : (Mi) 
where, 



(21) 
(22) 

(23) 



£±($) = (7*7*) « ± (7r7r)iJ ~ j ((e<et)« ± {e r e r )&) . 

and T± are formally given in Crittenden et al. (2002). The tangen- 
tial and radial ellipticity parameters et and e r are defined such that: 



et 



cos 2(f) sin 2(f) 
— sin 2(f) cos 2(f> 



ei 



(24) 



where (f> is defined to be the angle between the x axis and the line 
joining a galaxy pair. 

In the absence of contaminating non-lensing sources we ex- 
pect M± ~ 0, but this is found not to be the case in all the weak 
lensing surveys to date that have measured this diagnostic statistic. 
We therefore calculate the contribution to the aperture mass B mode 
from the four intrinsic alignment models defined in Sections 2.1 
and 2.2. This is calculated for two different depth weak lensing 
surveys, the RCS with a median redshift z m w 0.56, (Hoekstra 
et al. 2002), and the VIRMOS-DESCART survey with z m w 1.0, 
(Van Waerbeke et al. 2002). 

Fig. 5 shows the predicted angular shear correlation function, 



<77*>Wi 



|C7(i9), for each intrinsic alignment model cal- 



culated for both survey depths. We also plot for comparison the 
shear correlations expected from weak gravitational lensing, for an 
Q. m — 0.3, JIa = 0.7 cosmology with a T = 0.2, erg = 0.8, non- 
linear CDM matter power spectrum, Pg calculated using the fitting 
formula from Smith et al. (2003): 



(77*)(#)wgl 



1 

2^ 



dkkP K (k)J (M). 



(25) 



We calculate the aperture mass B mode contributions from the 
shear correlation functions shown in Fig. 5, through equation (22). 
CNPT find that {etet)$ ~ (e r e, )^ at small angular scales, deviat- 
ing slightly at separations i? > 3 arcmin, and this is consistent with 
the results we find. We therefore make the simplifying assumption 
that (etet)^ = (e r e r )^ at all angular scales, i.e that there is no pre- 
ferred intrinsic tangential or radial alignment of galaxy ellipticities, 
and hence = such that 



<Mi (IA) (0)> = 



d-dti 

A2 



(26) 



where £+(#) = \Ci{$) and we use the following analytic expres- 
sion for T+(x), derived by Schneider, Van Waerbeke & Mellier 
(2002), which vanishes at x > 2: 



T+(x) 



6(12- 15a; 2 ) 



2 x 
1 arcsin — 

7T 2 



+ 



cV4- 



IOOtt 

2\ 



x (110 + 2320x 2 - 754a; 4 + 132a; 6 - 9a; 2 ) . (27) 
Fig. 6 shows the aperture mass B mode measurements reported 
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Figure 5. Predicted intrinsic alignment angular shear correlation functions 
(77*) j, Jing (filled), HRH (stars), CNPT (circles) and HRH* (triangles). 
The intrinsic alignment contributions are compared to the angular shear cor- 
relations from ACDM weak gravitational lensing, normalised to erg = 0.8, 
(dashed), for the RCS with depth z m =0.56, (upper) and for the VIRMOS- 
DESCART survey with depth z m = 1.0, (lower). 



by the RCS (Hoekstra, Yee & Gladders 2002), (upper panel), and 
VIRMOS-DESCART survey, (Van Waerbeke et al. 2002), (lower 
panel), compared to the contributions from the four different intrin- 
sic alignment models. The observed B modes can be attributed in 
part to residual systematics remaining in the data and in part to in- 
trinsic galaxy alignments. Whilst most models for intrinsic galaxy 
alignments predict B modes there are alternatives such as the model 
proposed by Catelan, Kamionkowski & Blandford (2001) which 
does not. Known data systematics introduce positive correlations 
and we therefore consider these observations as upper limits for 
intrinsic alignment B modes within these surveys. This is proba- 
bly a valid assumption on small scales over which the distorting 
point spread function (PSF) is relatively constant, provided a con- 
sistent subtraction method for the PSF has been applied. Highly 
anisotropic PSF distortions which vary from e* — > — e* across 
the image could introduce a negative aperture mass B mode on 
large scales. This could in principle, introduce a systematic neg- 



RCS data 



1 




8 (arcmins) 

VIRMOS-DESCART Survey 

1 



1 












Jing 

- - HRH * 
-• - CNPT 
■■■■ HRH* - 














jsj { $f; 












) \ 


v 1 — . 
















■ ■ 1 


1 


1 1 L 


1 l_ 


1 1 1 







1 2 5 10 20 



8 (arcmins) 

Figure 6. Intrinsic alignment contributions to the aperture mass B mode 
{Ml); Jing(solid), HRH (dashed), CNPT(dot dash) and HRH* (dotted). 
The upper panel compares intrinsic alignment predictions for the RCS sur- 
vey with their measured B mode. The lower panel compares intrinsic align- 
ment predictions for the VIRMOS-DESCART with their measured B mode. 

ative aperture mass B mode at large scales > 10 arcmin, prevent- 
ing the use of B modes as an upper limit for intrinsic alignment 
contamination at these scales. 

With no evolution in galaxy clustering, considering the ob- 
served B modes as upper limits for 8 < 10 arcminutes, we find that 
the Jing and HRH models are strongly rejected by the RCS results 
and that the CNPT and HRH* models are favoured by both the 
VIRMOS-DESCART and the RCS. We note that none of the in- 
trinsic alignment models can account for the significant B mode 
at large angular scales 6 > 15 arcmin found in the VIRMOS- 
DESCART data, which may be caused by other data systematics 
(Hoekstra 2003). 

4 ESTIMATION AND REMOVAL OF THE INTRINSIC 
ALIGNMENT SIGNAL: SHEAR CORRELATIONS 

Decomposing cosmic shear correlation measurements into E and B 
modes provides an estimate of the intrinsic alignment contamina- 
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tion in weak lensing measurements. However it is not clear how to 
correct the E mode given a measured non-zero B mode. To date the 
practice has been to add in quadrature the B mode to the E mode 
error (Hoekstra, Yee & Gladders 2002). An alternative technique, 
applicable to weak lensing surveys with photometric redshift esti- 
mates, is to remove, essentially completely, the intrinsic alignment 
signal by downweighting pairs of galaxies which are likely to be 
physically close (King & Schneider 2002; Heymans & Heavens 
2003, hereafter HH). King & Schneider (2003) have also proposed 
separating cosmic shear from intrinsic galaxy alignments using cor- 
relation function tomography. 

In this Section, we extend the method of HH, to determine the 
intrinsic alignment signal, as well as to remove it from the shear 
correlation function. The method is reviewed in the appendix; close 
galaxy pairs are ignored in the shear correlation analysis, but in the 
optimised case, no more are removed than is strictly necessary, to 
avoid needlessly increasing the shot noise. 



4.1 Constraining intrinsic galaxy alignments with 
COMBO-17: Method 

COMBO- 17 is a deep multi-colour wide-field optical survey car- 
ried out with the Wide-field Imager (WFI) at the MPG/ESO 2.2m 
telescope that spans 1.25 square degrees to limiting magnitude 
R = 25.5 in five separate regions (Wolf et al. 2001). Multi-colour 
observations through a total of 17 filters, 5 broad-band (UBVRI) 
and 12 narrow band filters ranging from 420 to 914 nm, specifi- 
cally chosen to facilitate accurate photometric redshift estimation 
with errors A z < 0.05 reliable to R < 24 currently span 1 square 
degree of the survey area. The WFI instrument consists of a 4 x 2 
array of 2048 x 4096 CCDs with pixel scale 0.238" . Weak lensing 
studies have been carried out on the deep R-band images observed 
during the best seeing conditions (Brown et al. 2003), and it is a 
subset of this R-band selected galaxy sample, limited to R < 24 
with a median redshift z m ~ 0.6, that we will use in this analysis. 
The selected survey area totals 0.75 square degrees in the A901, 
CDFS and Sll fields, yielding a catalogue of 3.55 x 10 4 galaxies 
with KSB shape measurements and photometric redshifts accurate 
to A z = 0.042, (see Brown et al. 2003 for further details). 

The total observed shear correlation function in the COMBO- 
17 dataset can be expressed as 



C(9) = C lens (9) + d(9) + C aya (6), 



(28) 



where G\ C ns(9) is the cosmic shear weak lensing signal (77*} 
equation (25), and C\{9) is the signal due to the intrinsic alignment 
of galaxy shapes equation (5). C sy s(9) is the correlation function 
due to any observational systematic effects which may be present 
in the dataset (see Brown et al. 2003 for a robust estimate of sys- 
tematic errors in the data). C\{9) can effectively be eliminated by 
excluding galaxy pairs which are closer than a(9)A z in redshift, 
where A z = 0.042 is the typical redshift error for the COMBO- 
17 galaxies and a(9) is optimised to minimise the total error from 
intrinsic alignments and shot noise (see appendix). By measuring 
the shear correlation function for pairs of galaxies which are suf- 
ficiently distant from each other, we exclude intrinsic alignments. 
Furthermore, if we also measure the shear correlations for only the 
close galaxy pairs which have been excluded from the analysis, we 
have two independent correlation functions from which we can es- 
timate the intrinsic alignment signal from the close galaxy pairs: 

Cl A (9)= Cclose (0) - C dis tan t (9) - AC lcns (9), (29) 

where C c iose(0) is the correlation function measured from the close 



pairs, Cdistant(#) is the correlation function for the distant pairs 
and ACi ns(#) is the difference in the correlations caused by weak 
gravitational lensing for the distant and close pairs. Equation (29) 
assumes that in the distant data set the intrinsic alignment signal 
is zero and that C sys (9) is the same as for the close pairs. To en- 
sure all intrinsic alignment contamination has been removed from 
the distant dataset, we choose the most conservative estimate of 
7?jing(r) in our HH weighting scheme. The assumption for C ByB (9) 
is reasonable since there is no reason to believe that systematic ef- 
fects would depend on the redshift separation of galaxy pairs. To 
calculate ACi C ns(#), in Fig. 7 we plot the expected lensing sig- 
nal for a median redshift z m = 0.6 ACDM model, normalised 
to as — 0.8, for both the distant-pair dataset and for the close 
pairs, according to equation (A 13) in HH. Although the difference 
is small, at the 10 -5 level, and could be ignored, we will include it 
in our analysis assuming that ACi en s(#) is approximately the same 
for each field. The weak lensing signal Ci ens (0) will however dif- 
fer significantly between the three fields, A901, CDFS and Sll, 
due to the significant mass concentrations in two of the fields. The 
Sll field contains a fairly large cluster (Abell 1364) at a redshift 
of z — 0.11 while the A901 field includes a supercluster system 
(Abell 901/902) at z — 0.16. In the following analysis, we have 
therefore applied equation (29) to the three fields individually and 
combined the three resulting measurements of C\a(9) with mini- 
mum variance weighting, where the weight Wi for each measure- 
ment with associated error Oi is given by Wi = 1/of . This ensures 
that Ciens(#) is eliminated in the subtraction. Note also that, for the 
various C(9) calculations, we have removed cluster member galax- 
ies from the A901 and Sll fields to eliminate contamination of the 
field intrinsic alignment signal by intra-cluster galaxy alignments 
(Plionis et al. 2003). 

To estimate C c ioso (9) and Cdistant (9) for the three fields, we 
follow the analysis described in Brown et al. (2003). We split each 
field into eight chip-sized sections and calculate the relevant corre- 
lation function (ci jS ) for each section, 



<s,b(9) 



EabW»b7(fl»)7*(gfc) 



(30) 



where Wab is the weight for the galaxy pair ab. For the distant-pair 
data set i — 'distant', Wab = if \z a — Zb\ < aA z and W a b = 1 
otherwise. For the close pairs i — 'close' the weighting is simply 
reversed. Our estimate of the correlation function for the entire field 
is the average, 



(31) 



where N s is the number of sections in a field. Note that here we are 
making the approximation that sections of the same field are uncor- 
rected. Neighbouring sections of the same field are, in fact, corre- 
lated but our approximation will be valid on scales which are small 
compared to the size of the sections (= 8 x 16 arcmin). With this 
approximation, an estimate of the covariance matrix of the d(9) 
measurements is given by 



cov[a(0)Q(0')] ^ 



(32) 



We have investigated the sensitivity of the results to the choice of 
the number of sections each field is divided into, considering 4, 8 
and 16 sections. The most probable solution for a8(n m /0.27) ' 6 
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Figure 7. Predictions for the weak lensing shear correlation functions 
for the distant pairs and for the close pairs as described in the text. The 
difference in amplitude between the predicted signals is ACi ena (#) < 
1.5 X 10 — 5 at all measured scales. 



varies by about 0.03, much less than the statistical error, which is 
about 0. 1 (see Section 6). We do, however, add this error in quadra- 
ture in the quoted results. 

4.2 An observationally constrained intrinsic alignment model 

We have applied equations (29), (31) and (32) to the three fields, 
A901, CDFS and Sll and have combined the resulting measure- 
ments with minimum variance weighting to yield our final mea- 
surements of the intrinsic alignment signal for the close pairs in 
the COMBO-17 dataset. These final measurements are shown in 
Fig. 8 along with the predicted intrinsic alignment signal for the 
close pairs, 



Cia(0) 



J dz a dz h <t> z (z a )<t> z (z b ) [1 + igg] {W a b)r](r ab ) 

J dz a d,Z b (j) z (Za)(t>z{zb) [1 +€,gg] {W a b) 



(33) 



where we use three different models for the intrinsic alignment 
r](r ab ), and (W ab ) is given by 



(W ab ) = 



where, 



2^ 



dx e 



[erf (y + x) - erf (v + x)] (34) 



V2A Z V V2A Z V2A Z 

(Heymans & Heavens 2003). Fig. 9 shows the correlation matrix of 
the Cia{0) measurements, defined by 

cov(0,0') 



Zb - 



aA z 



Zb 



aA z 



(35) 



cor(f?, 



^cov(6>,6>)cov(6", 



(36) 



where cov((9, 9') — cov[Cia(#), Cia(0')] is the covariance matrix 
of the C\a{6) measurements, which we estimate from the data us- 
ing equations (32) and (29). One can immediately see from Fig. 8, 
in agreement with the aperture mass B mode analysis in Section 
3, that the measured COMBO-17 C\a(Q) strongly rejects the Jing 
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Figure 8. The intrinsic shear correlation functions from close galaxy pairs 
in the COMBO-17 survey, compared to predictions from Jing (solid), HRH 
(dashed) and HRH* (dotted). The CNPT prediction lies slightly above the 
HRH* prediction but on this scaling is indistinguishable and is therefore not 
plotted. 
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Figure 9. Correlation matrix of the optimally combined Cia(S) measure- 
ments plotted in Fig. 8. The area of each circle is proportional to the degree 
of correlation between points i and j . Filled circles denote that the points 
are correlated whereas unfilled circles denote an anti-correlation between 
the points. The points are the same as those plotted in Fig. 8, numbered 1 to 
10, in order from left to right. 



model for the intrinsic alignment signal while the HRH model is 
also highly inconsistent with the data. The modified HRH model, 
HRH*, introduced in Section 2.2, and therefore also the CNPT 
model, are much better fits to the data. We note also that the C\a{0) 
measurements are consistent with a null result on scales > 1 arcmin 
although, from Fig. 9 it's clear that there are some (anti-) correla- 
tions between the correlation function measurements at different 
scales. 

We calculate the likelihood of our HRH* model from the 
COMBO-17 results allowing the amplitude to vary, see Fig. 10. We 
find that the best fit is at zero amplitude, and that an HRH* model 



© 0000 RAS, MNRAS 000, 000-000 



Intrinsic Galaxy Alignments: estimation and removal 9 



-0.5 



1 - 



-1.5 



— CQMBO-17 

— COMBO- 17 (evo) 
SuperCOSMOS 

— HRH* 



0.002 



0.004 0.006 
Amplitude A 



0.008 



0.01 



Figure 10. The likelihood of the amplitude A of the model fit for the in- 
trinsic galaxy alignment correlation function equation (13), as determined 
from the COMBO-17 photometric redshift dataset assuming no evolution 
galaxy clustering, (solid), and assuming stable galaxy clustering, (dashed), 
see Section 5.2. Over-plotted is the likelihood of the amplitude from the 
modified HRH* results, (dot dash), and the likelihood of the amplitude from 
the SuperCOSMOS ellipticity variance measurements, (dotted), as detailed 
in Section 2.2. 



with an amplitude greater than A = 0.0054 can be rejected with 
95% confidence. We therefore find an observationally constrained 
maximum amplitude intrinsic alignment model such that: 

, , 0.0054 „„ „ , 
Vcn(r) < 5- 95yo confidence. (37) 



Note that using 16 x 16 arcminute sections in the correlation es- 
timator equation (31), we find consistent results with a 95% upper 
limit A = 0.0039. 



4.3 Application of the HH weighting scheme to COMBO-17 

With our constraints on the amplitude of the intrinsic alignment 
signal from the previous section, we can now apply the HH weight- 
ing scheme to the COMBO-17 data. In order to compare with the 
expected ACDM weak lensing shear signal, we apply the weight- 
ing scheme to the CDFS and S 1 1 fields only, in order to avoid the 
much larger shear signal seen in the A901 field due to the presence 
of the supercluster (Brown et al. 2003). To calculate the optimal 
a(8) weights, (see appendix), we use the COMBO-17 redshift dis- 
tribution for the CDFS and SI 1 fields, with our chosen input model 
for the intrinsic alignment signal as given by the upper limit of 
Vci7(r) equation (37), thereby ensuring removal of any feasible 
intrinsic alignment contamination. With values for a(8), shown in 
Fig. A2, we calculate the total and distant-pair correlation functions 
as given by equations (30) and (31), the results of which are shown 
in Fig. 11. The total C{9) (circles) will include intrinsic alignment 
contamination, whereas the distant-pair signal (triangles) excludes 
virtually all intrinsic alignment at the expense of a small increase 
in the shot noise. Over-plotted is the cosmic shear prediction for 
a ACDM model, equation (25), normalised to ag = 0.8, where 
the COMBO-17 photometric redshift distribution for the CDFS and 
S 1 1 fields has been used for this calculation. 
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Figure 11. Application of the intrinsic alignment suppression weighting 
scheme as described in the text. The circles are the total correlation func- 
tion C{9) measurements from all galaxy pairs whereas the triangles are 
the C(9) measurements from just the distant galaxy pairs. The curve plot- 
ted is the expected weak lensing signal for a ACDM model normalised to 
a s = 0.8. 



This result is a good example of the application of the HH 
weighting to a weak lensing data set, effectively removing the sys- 
tematic error from intrinsic galaxy alignments whilst producing 
only a small increase in the shot noise. Due to the fact that we are 
only using a small subset of the COMBO-17 data set (0.5 square de- 
grees limited to magnitude R < 24), we are unable to use this noisy 
result directly to investigate the effect on cosmological parameter 
estimation. We can, however, use the likelihood of the amplitude 
of the intrinsic alignment signal to calculate its effect on parameter 
estimation from the full COMBO-17 deep i?-band sample. We do 
this in Section 6. 



5 EFFECT OF GALAXY CLUSTERING EVOLUTION 

All previous studies of intrinsic galaxy alignment have assumed 
weak evolution in galaxy clustering which becomes important 
when converting measured 3D intrinsic ellipticity correlations, 
7/(r), into angular correlation functions Ci(9), see equation (5). 
If there is evolution in galaxy clustering such that at high redshift 
there are proportionally less nearby galaxy pairs than at low red- 
shifts, this will lower predictions of Ci{9). The clustering evolu- 
tion of galaxies has often been quantified, although without strong 
theoretical motivation, by a redshift-dependent two point correla- 
tion function of the form (Groth & Peebles 1997) 



t M {r,z)= (-Y(l + zy 
\roJ 



(3+e) 



(38) 



Observational studies have found that this simple redshift- 
dependent model provides a poor fit to data (Wilson 2003; Mc- 
Cracken et al. 2001), which is most likely due to the fact that differ- 
ent populations of galaxies have different intrinsic clustering. Pa- 
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rameters for equation (38) are therefore observationally challeng- 
ing to determine and as yet fairly uncertain. Wilson (2003) investi- 
gates the redshift evolution for a single population of galaxies in the 
UH8K weak lensing fields and it is these results that we will use in 
order to investigate the effects of galaxy clustering evolution in in- 
trinsic alignment studies such that 7 = —1.8, ro = 5.25/i _1 Mpc 
and e = 0. For e — 0, the evolution is of stable clustering, where 
the galaxies are dynamically bound and stable at small scales. 

Fig. 12 shows the predicted angular shear correlation function, 
(77*} 1 — \Ci{9) calculated including galaxy clustering evolu- 
tion for each intrinsic alignment model for the RCS and VIRMOS- 
DESCART surveys. We also plot for comparison the shear cor- 
relations expected from weak gravitational lensing, equation (25). 
Comparing this with Fig. 5 we note that by assuming stable cluster- 
ing evolution, the intrinsic alignment contribution at small angular 
scales is significantly reduced compared to previous results. This 
is due to the fact that it is at small angular scales that high redshift 
galaxy pairs are close enough to contribute to any intrinsic corre- 
lations. If galaxies are less strongly clustered at high redshift then 
there will be fewer high-redshift close pairs and hence the angular 
intrinsic correlation function is reduced. 



5.1 Aperture mass B modes 

We repeat the analysis of Section 3 with the intrinsic alignment cor- 
relation functions now calculated including galaxy clustering evo- 
lution, as shown in Fig. 12. Fig. 13 shows these results, comparing 
the aperture mass B mode measurements reported by the RCS, (up- 
per panel), and VIRMOS-DESCART survey, (lower panel), and the 
four different intrinsic alignment B modes. Unlike the clear result 
from Section 3 the inclusion of stable galaxy clustering reduces the 
expected B mode contribution from intrinsic galaxy alignments, 
permitting all intrinsic alignment models where the observed B 
modes are considered as upper limits. 



5.2 COMBO-17 Correlation Analysis 

We have also investigated the effect of including galaxy evolution 
in the intrinsic alignment models for the COMBO-17 correlation 
analysis. The resulting Cia(#) measurements are shown in Fig. 14 
along with the predictions for the Jing, HRH and HRH* intrin- 
sic alignment models where stable clustering galaxy evolution has 
been included. Here, we find that, although not as clear a result as 
found in Section 4, we can still exclude the Jing model and, once 
again, the HRH* model is favoured. We calculate the likelihood 
of our HRH* model which includes clustering evolution, allowing 
the amplitude A to vary, see Fig. 10. We find again that the best 
fit is at zero amplitude, and that an HRH* evolving model with 
an amplitude greater than A = 0.0096 can be rejected with 95% 
confidence. 



6 IMPLICATIONS FOR WEAK LENSING 
MEASUREMENTS 

6.1 Cosmological parameter constraints 

We now turn to the implications of our intrinsic alignment am- 
plitude constraints for cosmological parameter estimation. To do 
this we follow Brown et al. (2003) using measurements from the 
deep R < 25.5 COMBO-17 data of the two correlation functions, 

Ci(6>) = £[7t7t : 6>] and C 2 (0) = E[-y r -y r : 8} from equation 
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Figure 12. Predicted intrinsic alignment angular shear correlation functions 
(77*)/ including stable galaxy clustering evolution, Jing (filled), HRH 
(stars), CNPT (circles) and HRH* (triangles). The intrinsic alignment con- 
tributions are compared to the angular shear correlations from ACDM weak 
gravitational lensing, normalised to erg, = 0.8, (dashed), for the RCS with 
depth z m = 0.56, (upper) and for the VIRMOS-DESCART survey with 
depth z m = 1.0, (lower). 



(31), with W a b = 1 for all galaxy pairs, in order to obtain a joint 
measurement of the normalisation of the matter power spectrum, 
as and the matter density, Q m . We perform a x 2 fitting procedure 
on the correlation function measurements ordered in a data vector 
d = {Ci(e 1 ),...,Ci(6 n ),C2{ei),...,C 2 {e„)}, for a set of theo- 
retical parameters, (as, f2 m , Ho, A), where Ho is the Hubble pa- 
rameter and A is the amplitude of the intrinsic alignment signal in 
equation (13), by calculating 



X 



[d-x(a s ,n m ,H ,A)] T V- 1 [d-x(as,n m ,H ,A)],(39) 



where x(crg, Q m , Ho, A) is a theory vector containing the C\{8) 
and C2(d) correlation functions calculated for the cosmological 
model. V = (dd T ) is the covariance matrix of the data measure- 
ments which can be estimated, assuming that sections of the same 
field are uncorrelated, from equation 32. We assume a flat cosmol- 
ogy Q m + Qa ~ 1 and a nonlinear CDM matter power spectrum, 
Ps, calculated using the fitting formula from Smith et al. (2003) 



© 0000 RAS, MNRAS 000, 000-000 



Intrinsic Galaxy Alignments: estimation and removal 1 1 



RCS data 



<u 

T3 
O 

E 

CO o 
CO 

o 
E 

t U 

a) o 

< 



Jing 

- - HRH 

- ■ - CNPT 
.... HRH* 



■ ■ ■ i 



1 2 5 10 20 

6 (arcmins) 



I 

o 



I 1 1 1 1 1 — I — I — [- 



— Jing 

- - HRH 
••HRH* 



2 5 
6 (arcmins) 



10 



VIRMOS-DESCART Survey 



E 

DD 



E 

a) o 

< 



: 1 1 










Jing 












- - HRH : 












- ■ - CNPT . 












•■■• HRH* ; 






































:■;■!' 




■ ■ ■ 


■ ■ 


i i i 





2 5 10 20 

6 (arcmins) 



Figure 13. Intrinsic alignment contributions to the aperture mass B mode 
{Mj) assuming evolution in galaxy clustering; Jing(solid), HRH (dashed) 
and CNPT(dot dash). The upper panel compares stable clustering intrinsic 
alignment predictions for the RCS survey with their measured B mode. The 
lower panel compares stable clustering intrinsic alignment predictions for 
the VIRMOS-DESCART with their measured B mode. 



with the transfer function of Eisenstein & Hu (1999), where the 
initial slope of the power spectrum is n = 1. In contrast to (Brown 
et al. 2003), we do not fix the Hubble constant Ho, but marginalise 
over it, with a Gaussian prior p(Ho) set by the WMAP results with 
H = 72 ± 5 kms^ 1 MpcT 1 (Spergel et al. 2003). 

We calculate \ 2 equation (39), for values of cr 8 ranging from 
0.3 to 1.5 and values of Q. m ranging from 0.1 to 1.0, on a grid with 
0.05 spacing in both these parameters, for two cases (a) the case 
where we ignore any effect of intrinsic alignments setting A = 
and (b) the case where we also vary values of A from 0.0 to 0.01 in 
steps of 0.0005 and then marginalise over the intrinsic alignment 
amplitude. For the marginalisation we use the probability distri- 
bution for the intrinsic alignment amplitude p(A) obtained from 
COMBO- 17 and shown in Fig. 10 such that 

p(a 8 ,Q m )= I dA I dHop(as,n m ,H ,A)p(Ho)p(A), (40) 



Figure 14. The intrinsic shear correlation function from close galaxy pairs 
in the COMBO- 17 survey, compared to predictions from Jing (solid), HRH 
(dashed) and HRH* (dotted) for stable clustering galaxy evolution. 



where p(os, fi m , Ho, A) oc exp(— x 2 /2). As in Section 3, we as- 
sume that intrinsic galaxy alignments have no preferred tangential 
or radial alignment, and that d(8)u ~ Ci(9) rr w ^Ci(0), where 
Ci(6») is now calculated for the deep R < 25.5 COMBO-17 esti- 
mated redshift distribution detailed in Brown et al. (2003). 

The resulting constraints in the <j 8 — f2 m plane are shown in 
Fig. 15. Assuming no intrinsic alignment signal, we obtain a best 
fit measurement for the normalisation of the power spectrum of 



a 8 (n m /0.27) = 0-74 ± 0.11, 



(41) 



while for the case where we have marginalised over the intrinsic 
alignment amplitude derived assuming no clustering evolution, we 
find our estimate drops to 



a 8 (fi m /0.27) ' 6 = 0.71 ± 0.11. 



(42) 



The errors are larger than quoted in Brown et al. (2003) be- 
cause we have marginalised over Ho rather than fixing it at 68 
km s _1 Mpc _1 , but they remain relatively small in comparison to 
other surveys, bearing in mind the COMBO-17 survey area. This 
is a result of the good determination of source redshifts in the 
COMBO-17 survey, and the relatively high number density of re- 
solved background sources used in this analysis. If we marginalise 
over the Ho distribution obtained from CMB and 2dF with Ho = 
68 ± 5 km s" 1 Mpc" 1 (Efstathiou et al. 2002; Percival et al. 2002), 
these numbers change to a 8 (r2 m /0.27) ' 6 = 0.78 and 0.77, with 
an error of 0.12. Note that all these error bars include a small error 
added in quadrature from the imperfect estimation of the covari- 
ance matrix (section 4.1). 

Note that if we instead use measurements of the total corre- 
lation function, C{6) = Ci(6) + C2(9) we find the same reduc- 
tion of 0.03 in our estimates of <r 8 . With the intrinsic alignment 
constraints from the SuperCOSMOS survey we find a reduction in 
our estimates of crs of 0.01. Fig. 16 shows the COMBO-17 mea- 
surement of the shear correlation function along with the predicted 
lensing shear correlation function with the HRH* intrinsic align- 
ment signal subtracted. 

For the case of stable clustering galaxy evolution, the resulting 
angular intrinsic alignment contamination significantly decreases, 
and we find that the constraints in the as — Q. m plane are unchanged 
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model j7hrh* (r ) we find that 

4^y ( g= larcmin) = 0.02 (CFHTLS), 
= 0.01 (SNAP), 



(43) 



Figure 15. The probability surface for erg and f2 m from COMBO-17 as 
calculated using the shear correlation function measurements C'i (9) and 
C'2(8) in combination, after marginalising over the Hubble constant. The 
intrinsic alignment signal is either assumed to be zero (solid contours) or is 
marginalised over (dashed contours). The inner and outer contours in each 
case correspond to A\ 2 = 2.3 and 6.17. confidence regions respectively. 



n =0.3, ct„ = 0.75, z =0.85 

m o m 




6 (arcmin) 

Figure 16. The total shear correlation function for the deep COMBO-17 
survey, along with the best fit ACDM weak lensing correlation function as- 
suming no intrinsic alignments (solid). Subtracting off the the best-fitting 
intrinsic alignment model (HRH*) decreases the lensing correlation func- 
tion to the dashed line, and with the 95% upper limit (obtained from the 
shallower survey with photometric redshifts) the lensing correlation func- 
tion decreases to the dotted line. 



after marginalising over the intrinsic alignment signal, derived as- 
suming stable clustering galaxy evolution. 



6.2 Implications for future surveys : CFHTLS and SNAP 

The implications for weak lensing surveys is obviously dependent, 
as shown in Section 5, on our understanding of galaxy cluster- 
ing evolution. Until agreement is reached upon the true redshift 
dependence of galaxy clustering we consider the contamination 
to weak lensing measurements derived assuming no evolution in 
galaxy clustering, which should be considered as upper limits as 
redshift evolution will dilute the signal. For the modified HRH 



where we have calculated Ci(6) from equation (5) with <f>z(z) as 
given in equation (6), with z m = 1.0 for the CFHTLS wide survey, 
and z m — 1.23 for the SNAP wide survey. (77*) has been calcu- 
lated from equation (25) with £7 m = 0.3, Qa = 0.7, as = 0.8 and 
T = 0.2. The contamination of weak lensing measurements at an- 
gular scales 6 > 1 arcmin is therefore less than 2% for CFHTLS, 
and less than 1% for SNAP. These estimates are essentially un- 
changed if we use the SuperCOSMOS measurement of intrinsic 
alignments. Using the 95% upper limit for the intrinsic alignment 
amplitude r]cn(r), obtained from the COMBO-17 data, these es- 
timates are multiplied by approximately 5, giving the 95% upper 
limits on any possible intrinsic alignment contamination at 10% for 
CFHTLS and 5% for SNAP For comparison, with stable clustering 
evolution, contamination for both surveys is practically negligible, 
less that 0.2% for CFHTLS, and less that 0.1% for SNAP. 

Refregier et al. (2003) have shown that for SNAP, the appli- 
cation of redshift tomography as proposed by Hu (1999), increases 
the accuracy of cosmological parameter estimates by a factor of 
about 2. This method involves splitting the galaxy sample into two 
or three redshift slices, and is therefore more susceptible to higher 
levels of intrinsic alignment contamination (Croft & Metzler 2000). 
Using the redshift slice distributions as proposed by Refregier et al. 
(2003), we have calculated the intrinsic alignment contamination 
for the best-fit HRH* model, and found that for the three SNAP red- 
shift slices with median redshifts 0.8, 1.3 and 1.9, the contamina- 
tion at angular scales 6 > 1 arcmin is found to be, 7.0%, 2.5% and 
0.7% respectively. For the lowest redshift bin, this level of contam- 
ination is significant, especially if we conservatively consider the 
95% upper limit from our COMBO-17 constraint, which increases 
the contamination by a factor of approximately 5 to ~ 35%. It 
will therefore be important, when using tomography in weak lens- 
ing analysis, to apply a close galaxy pair downweighting scheme in 
order to remove this systematic error. Note that when considering 
stable galaxy clustering evolution the contamination is 1% for the 
lowest SNAP redshift bin. 



7 CONCLUSIONS 

The weak correlation of the ellipticities of galaxy images is an in- 
dicator of gravitational lensing and a powerful tool for the study of 
dark matter on large scales. In this paper we have used two meth- 
ods to estimate the extent to which this statistic is contaminated 
by the intrinsic physical alignments of galaxies. Our main conclu- 
sion is that the effect is relatively small, but not entirely negligible, 
and, for the COMBO-17 survey, leads to a reduction in the derived 
amplitude of mass clustering of around 3%. 

For the brighter part (R < 24) of the COMBO-17 survey, 
which has photometric redshifts, we removed close pairs from the 
shear correlation analysis, removing the intrinsic alignment sig- 
nal, as described by Heymans & Heavens (2003). Comparing this 
shear correlation function from the distant pairs with that of the 
close pairs then allows us to estimate the intrinsic alignment signal. 
We have also placed limits on the intrinsic alignment signal from 
analysis of the aperture mass B mode in the RCS and VIRMOS- 
DESCART surveys. We find a consistent picture that this signal is 
lower than that expected from analysis of numerical simulations 
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(Heavens, Refregier & Heymans 2000; Jing 2002), but in broad 
agreement with the semi-analytic calculation of Crittenden et al. 
(2002). 

We have reanalysed the numerical simulations to include two 
effects which were originally ignored. These are a misalignment 
between the angular momentum of the baryons and the halo (van 
den Bosch et al. 2002), and the finite thickness of disk galax- 
ies (Crittenden et al. 2002). Both these effects reduce the intrin- 
sic galaxy ellipticity correlation function to a level similar to that 
found by Crittenden et al. (2002), and consistent with the level de- 
termined observationally in this paper and the level measured in the 
SuperCOSMOS survey (Brown et al. 2002). Note that other effects 
such as gas-dynamical interaction have not been taken into account, 
which could cause the intrinsic alignment signal to be lowered still 
further. 

Having estimated the contribution of intrinsic alignments to 
the brighter part of the COMBO-17 data, we have computed the 
likelihood for the contamination of the whole COMBO-17 sample 
which extends to R < 25.5. We compute the probability distribu- 
tion for <T8(fi m /0.27) ' 6 , now marginalising over the Hubble con- 
stant. Marginalising over the amplitude of the intrinsic alignment 
contamination, the mass clustering amplitude O8(£2m/0.27) ' 6 = 
0.71 ± 0.11. Ignoring intrinsic alignments leads to a systematic 
overestimate by 0.03. 

From the COMBO-17 results we have also calculated 95% 
upper limits for the expected contamination of future surveys, 
CFHTLS « 10.0%) and SNAP « 5%). With our theoretical 
model, or with the intrinsic alignment signal estimated from the 
SuperCOSMOS data, these predicted limits become 2% and 1% 
respectively. With current surveys these levels of contamination 
are small enough to be neglected, but will be significant in the 
error budget of future high-precision weak lensing surveys. Both 
CFHTLS and SNAP aim to produce accurate photometric redshift 
estimates for their galaxy sample enabling the use of redshift to- 
mography to further improve cosmological parameter estimation. 
This technique is susceptible to significant contamination from in- 
trinsic galaxy alignments due to the thin widths of the redshift 
bins, increasing the proportion of nearby galaxy pairs. We have 
shown that, even with the low amplitude intrinsic alignment model 
J?hrh*(V), with the proposed redshift distributions for the SNAP 
tomographic redshift bins, the lowest redshift bin will suffer con- 
tamination ~ 7%, with a 95% upper limit of 35% if we consider 
our observationally constrained upper limits for r)ci7{r). Since 
these surveys will have photometric redshift information, it will 
therefore be vital to remove the intrinsic alignment signal using 
the exclusion of nearby galaxy pairs as proposed by Heymans & 
Heavens (2003) and King & Schneider (2002). 

Our conclusions are affected by the clustering strength of 
galaxies, as this partly determines how many pairs of galaxies 
which are close on the sky are actually physically close together, 
and susceptible to physical interactions which could lead to intrin- 
sic alignments. The results we have presented so far assume that 
clustering is independent of redshift, but for illustration we have 
also investigated, without strong theoretical motivation, an evolu- 
tionary model corresponding to stable galaxy clustering. In this 
case, the effects at high redshift are reduced to a negligible level 
for COMBO-17, CFHTLS and SNAP, but could still be important 
in the case of weak lensing tomography analysis. 

In the process of applying weighting schemes as proposed by 
Heymans & Heavens (2003) and King & Schneider (2002) to future 
weak lensing surveys, there will potentially be some very interest- 
ing by-products. For example, with large area redshift slices, the 



method detailed in Section 4.1 could be applied in order to deter- 
mine the strengths of intrinsic galaxy alignments as a function of 
redshift, which throughout this paper we have assumed to be con- 
stant. In principle this could be a useful constraint for galaxy for- 
mation and evolution studies. With the SuperCOSMOS results and 
our COMBO-17 intrinsic alignment constraint favouring an intrin- 
sic alignment model which includes misalignments between baryon 
and halo angular momentum, there is now observational evidence 
indirectly supporting the finding by van den Bosch et al. (2002), 
which has important implications for formation of disk galaxies. 
We therefore conclude that although, for studies of weak gravi- 
tational lensing, the presence of intrinsic galaxy alignments is an 
inconvenience, they are an interesting topic in their own right. 
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APPENDIX A: OPTIMAL GALAXY PAIR WEIGHTING 

For a pair of galaxies with estimated redshifts z a and zt and associ- 
ated errors A z , following Heymans & Heavens (2003), we assign a 
zero weight if \z a — &b\ < ctA z and a weight of one otherwise. We 
choose q to minimise the total error on the shear correlation func- 
tion which has two components, one systematic error from intrin- 
sic galaxy alignments, oia, and one random error from shot noise. 
The optimum a value will depend on angular separation 9, the 
COMBO-17 redshift distribution, <f> z (z), shown in Fig. Al, and the 
median photometric redshift accuracy, A z — 0.042. The shot noise 
comes from the COMBO-17 intrinsic distribution of galaxy ellip- 
ticities, a e — 0.67. The estimated value for ctia depends on which 
r/(r) model we choose. The conservative approach is to consider 
the highest amplitude intrinsic alignment model, thereby ensuring 
all possible pairs of galaxies that could contribute to the contam- 
inating intrinsic alignment signal are removed. It is this approach 
that we use in Section 4. 1 in order to observationally constrain the 
intrinsic alignment contribution to COMBO-17. By doing this how- 
ever there is the expense of a potentially needless high residual shot 
noise, and therefore the preferred approach which we use in Sec- 
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Figure Al. COMBO-17 total redshift distribution for the A901, CDFS and 
S 1 1 fields, limited to R < 24, with known galaxy cluster members re- 
moved. 

tion 4.3, is to use the upper limit of r]cn(r) in order to optimally 
remove the intrinsic alignment contamination from the weak lens- 
ing signal in COMBO-17. 

Fig. A2 shows the optimal a values calculated for the 
COMBO-17 redshift distribution for the highest amplitude intrin- 
sic alignment model: Jing (circles), and the upper limit COMBO- 
17 model (stars). Spikes in the COMBO-17 redshift distribution, 
shown in Fig. Al, mean that a(0) is not a smooth function. The 
higher amplitude Jing intrinsic alignment model requires more 
galaxy pairs to be rejected than for the lower amplitude COMBO- 
17 model, where the intrinsic alignment contribution is of similar 
amplitude to the shot noise. In this case the best results are ob- 
tained with the rejection of only the closest galaxy pairs but this 
would not necessarily be the case for large area surveys where high 
galaxy number counts reduce the total shot noise. 

We derive a values for the Jing model considering no evo- 
lution in galaxy clustering (solid) for Section 4.2 and considering 
stable clustering (dashed) for Section 5.2. The angular correlation 
signal from intrinsic galaxy alignments are less when we include 
galaxy clustering evolution, therefore the optimal a values are less. 
The expected intrinsic alignment ellipticity correlation signal with 
stable clustering, is fairly constant for angular scales 8 < 10 ar- 
cmin and this is reflected in the optimal a values. As the shot noise 
decreases with increasing 8, the optimal a can increase to remove 
more of the intrinsic alignment signal without increasing the total 
error. 



Optimal a values 
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Figure A2. Optimal a values for the COMBO-17 redshift distribution for 
the Jing intrinsic alignment model (circles), and the upper limit from the 
COMBO-17 constrained model (stars). For the Jing model we show the 
optimal a values considering two different galaxy clustering models: zero 
evolution (solid) and stable clustering (dashed) 
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